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Abstract 

We investigate the spectrum of the free Neuberger-Dirac operator Z?n on the 
Schrodinger functional (SF). We check that the lowest few eigen- values of the Her- 
mitian operator D^D^ in unit of L~ 2 converge to the continuum limit properly. 
ON ! We also perform a one-loop calculation of the SF coupling, and then check the uni- 



versality and investigate lattice artifacts of the step scaling function. It turns out 
that the lattice artifacts for the Neuberger-Dirac operator are comparable in those 
of the clover action. 



1 Introduction 



Chiral symmetry plays an important role in the understanding of the strong inter- 
action. A solution to a realization of the exact chiral symmetry on the lattice is 
proposed by Neuberger [I] as the Neuberger-Dirac operator (overlap). Recently, 
dynamical overlap lattice QCD simulations started in Ref. [2] and see [3] for an 
overview of recent progress. Nowadays, thanks to developments of algorithms and 
powerful current computers, large scale simulations are feasible as shown by the 
JLQCD collaboration in Ref. [4J. In that course, even after removing systematic 
errors 0, finite size effects, cutoff effects and an ambiguity of chiral extrapolation, 
non-perturbative renormalization becomes an essential element for accurate quan- 
titative predictions. 

One elegant solution to this issue is the Schrodinger functional (SF) scheme [5] 
which is an intermediate scheme connecting the perturbative and hadronic regime. 
This method was shown to be useful to study the non-perturbative evolution over a 
wide range for various quantities, the coupling constant [6], the quark masses [7], the 
structure function [8j [9] , and the weak matrix elements [10] . By making use of the 



1 Of course, quenching is also one of the main sources of systematic errors. 
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scaling technique, one can complete a perturbative matching safely at relatively high 
energy, and then the renormalization group invariant quantities, like the lambda 
parameter [llj . the masses for the light quark [J2] and the heavy quark by HQET 
|13j . are determined without worrying about the systematic error of truncation of 
the perturbative expansion. Although so far the method was mainly used for Wilson 
type fermions [T4], there are several attempts for other fermion formulations, like 
staggered fermions |15l [TBI \T7\ and domain wall fermions [18] . Recently, formalisms 
for the Neuberger-Dirac operator on the SF have been proposed by Taniguchi [19] 
and Liischer [20J. The former employs the orbifolding technique and the latter 
is based on universality considerations. In Ref. [21[ I22j . the latter formulation is 
examined in the framework of the Gross-Neveu model. Here in this paper, we 
implement the formulation in QCD and study the spectrum of the free operator 
and calculate the SF coupling to one-loop order. Furthermore we investigate lattice 
artifacts of the step scaling function by comparing with the clover fermion action. 

The rest of the paper is organized as follows. In Section 12.11 we summarize the 
definition of the Neuberger-Dirac operator on the SF which is given by Liischer, and 
give some practical details about building the operator. Then we show the spectrum 
of the free Neuberger-Dirac operator in Section [31 In Section HI we present results 
for the fermion part of the SF coupling to one-loop order by making use of the 
Neuberger-Dirac operator. Furthermore we investigate lattice artifacts of the step 
scaling function in Section Finally, we give some concluding remarks and outlook 
in Section [6j In Appendices, we give an explicit form of the Wilson-Dirac operator 
on the SF in a time- momentum space for later use (}A"j) . and a discussion about a 
determination of a boundary coefficient at tree level ([B]), and we summarize some 
tables of numerical results (ICl) , 

2 Neuberger-Dirac operator on the SF 
2.1 Definition 

A massless Neuberger-Dirac operator on the SF (with size T x L 3 ) [20] is defined 
by 

= \{l- l -(U + U)}, (1) 

U = AX-^ 2 ,X = A^A + caP, (2) 
U = 75^75, (3) 

with a = a/(l + s). The operator follows the modified Ginsberg- Wilson (GW) 
relation 

7s£>N + £>n75 = aD N -/ 5 D N + A B , (4) 

where Ab is the exponentially local operator. The A in the kernel operator X of 
the inverse square root is given by 

A = l + s-aD w , (5) 

where D w is the massless Wilson fermion on the SF [14J. The tunable parameter 
s is taken in a range —0.6 < s < 0.6 in the following. The boundary coefficient c 



2 



represents the strength of the boundary operator P which is given by 



P 



-S Xj yS XO) y {S XOta P-. + S X0) T— a P+}- 



(6) 



In Ref. [20], it is shown that the kernel operator X = A^A + caP is bounded 
from below by the spectral gap of A'A on the infinite lattice if c > 1 holds, and 
furthermore it is mentioned that 



1 + s, 



(7) 



is the nearly optimal choice in order to achieve tree level 0(a) improvement. We 
investigate this point in some detail in Appendix [Bj and we conclude that, for the 
precision of our calculation here (and maybe for future simulations), this formula is 
accurate enough. Therefore in the following calculations, apart from the appendix, 
we alway set c to the value in eq.(|7|). In this paper, we restrict ourself to the massless 
case of the Neuberger-Dirac operator. 

To carry out perturbation calculations, it usually is convenient to move to mo- 
mentum space. In the SF setup, however, the Fourier transformation can be done 
only in the spatial directions, due to a lack of translation invariance for the time 
direction. Therefore, we work in a time-momentum space, 



a 2 e 

X 



-ipx 



(8) 



The explicit expression of the free D w in time-momentum space with the Abelian 
background gauge field in the SU(3) group [6j is shown in AppendixEJ Here we show 
a matrix expression of A in time-momentum space for a fixed spatial momentum p 
and a color b(= 1,2,3), 



g b (p;a) P_ 
P+ 5 6 (p;2a) 
























p+ 






5 6 (p;3a) P_ 

P+ '■■ '■■ 

P + 




g b (p;T-2a) 










g b (p;T-a) 



(9) 

which is a totally 4(T/a — 1) x 4(T/a — 1) matrix. Block elements P± = (1 ± 7o)/2 
and g b (p,xo) have Dirac spinor structure, and the latter is given by 



^3 3 
g b (p;x ) = s - -^2q b k (x ) 2 - i^2q b k (x )j k 



(10) 



k=l 



k=l 



where q and q are function of the spatial momentum p, the time xo and 9, and 
they are defined in Appendix [A] The dependence of b is caused by the color diag- 
onal background field. The angle 9 comes from the generalized periodic boundary 
condition for the spatial directions, 



ip(x + Lk) = e ie ij){x), ip(x + Lk) 



-i0 



ij){x) 



(11) 
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for k = 1,2,3. In time-momentum space, the boundary operator in eq.( 
sented as ^ 

Ptm = -5x ,y {$x ,aP~ + 5x ,T-a,P+}, 

and it's matrix expression is given by 

P_ 

Ptm 



is repre- 
(12) 












P + 



(13) 



As it is clear from the explicit form, even in the free case we cannot have an 
analytic form of the Neuberger-Dirac operator on the SF, due to the presence of 
the background field. Therefore we have to rely on an approximation to the inverse 
square root, even for perturbative calculations. We will return to this issue how to 
build the operator in Section | 



2.2 Distribution of e 

When one approximates the X -1 / 2 by a polynomial of X following Ref. [23) . one 
needs information about the lower u and upper v bound of X. When the ratio u/v 
is not too small, one can obtain its approximation with lower degree. Usually, u 
and v are set to the values of minimal and maximal eigen- value (or norm) of the X, 
therefore it is important to know the spectrum of X. Here we discuss the spectrum 
of X at tree level in the presence of the background field. 

At fixed momentum p and color b, we evaluate the minimal eigen-value and the 
norm of the kernel operator 

X\p) = (A\p)^A b (p)+caP tm , (14) 

and set them as the lower and upper range of the approximation [u b (p) , v b (p)] 

u b (p) = A min (X b (p)), (15) 
v b (p) = \\X b (p)\\. (16) 

Essentially, e b (p) = u b (p)/v b (p) controls the cost of the computation, since for a 
given precision it determines the degree (N) of the approximation polynomial. 

Its distribution over p and color b in the case of = with the non-zero back- 
ground field are shown in Figure Q] for several lattice sizes L/a = T/a = 6,12,24 
and s parameters s = 0.0, 0.5. As you can see in the figure, e is distributed in a rela- 
tively large range 0.01 < e < 1, therefore we concluded that it is better to calculate 
coefficients of the polynomial expansion for each p and b for saving time. The cal- 
culation of the coefficients is much cheaper than summing up the polynomial expan- 
sion. (You can use common lower and upper bound u comm = min p {, {A min (X 6 (p))} 
and Wcomm = max Pi 6 ||X 6 (p)|| for all momentum and color sector, but this is clearly 
inefficient.) 
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Figure 1: We show the distribution of e 6 (p) = u b (p)/v b (p) over the all spatial momenta 
p and the color b = 1, 2, 3 for the case of 9 = with the non-zero background field. We 
show here for lattice sizes L/a = T/a = 6, 12, 24 (from top to bottom) and s = 0.0 (Left 
panels) and s = 0.5 (Right panels). Due to the permutation symmetry for the spatial 
momentum, we only have to investigate a momentum configurations p% < P2 < Pi- We 
did not correct the weight in the histogram, therefore the distribution shown here is not 
fully correct one. However the distribution is practically meaningful, since we do not do 
a computation on a full configurations with size (L/a) 3 , but only on the configurations 
Ps < P2 < Pi (about (L/a f/6). 

2.3 How to build the Neuberger-Dirac operator 

To build the inverse square root in the Neuberger-Dirac operator, we adopt the 
Chebyshev polynomial approximation, 

N 

f(X) = X- 1 / 2 ~ f N (X) = Y / c k T k (Y), (17) 

fc=0 
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where Y = [2X — (v + u)]/(v — u) for the lower u and upper v bound of the X, and 
7*. is the Chebyshev polynomial of degree k. Main tasks to obtain the operator are 
two-fold: the first is a computation of the coefficients and the second is to sum 
up in eq,(|17p. Concerning the latter part, we use the Clenshaw summation scheme 
to maintain precision. For the former part, we examine two methods to compute 
the Cfc, the Remez algorithm and the Chebyshev interpolation, in order to check 
rounding off errors occurring in the perturbative calculation in Section [H 

Following Ref . [23] , we implement the Remez algorithm to obtain minimax poly- 
nomials. On the other hand, we follow the Numerical Recipes for the Chebyshev 
interpolation. In calculations of the spectrum of the Neuberger-Dirac operator in 
Section [3j we only use the Chebyshev polynomial method. While for computations 
of the one- loop coefficient of the SF coupling in Section [3J we employ both methods 
and show the results to digits which are common in both methods. 

When approximating the inverse square root of X, we demand a precision 



max 

u<x<v 



f(x) - /at (x) 



< 1(T 13 , (18) 



m 

and we alway check consistency when it is built 

X x (X- 1 / 2 ) 2 = /, (19) 
and we observe that errors on the right-hand side are less than 2 x 10 -13 . 



3 Spectrum of the free Neuberger-Dirac op- 
erator 

3.1 Spectrum of Z) N 

Since the U in eq.([2]) is not a unitary matrix, there is no guarantee that its spectrum 
is distributed on a unit circle whose origin is (1,0) as in the case of the infinite 
volume. However on the SF, as it is shown in Ref. [20], since 

H^ll = H^ll < 1, (20) 

one can see that 

[|oDn-1|| = -\\U + U\\ < 1. (21) 

Thus the spectrum of aD^ is contained in a unit disk which is enclosed by the unit 
circle. 

We show the actual distribution of the spectrum in the free case in Figure [2] 
for s = 0.0 and s = 0.5. In the figures, we only show 9 = case, but 9 = ir/5 
is also computed and has a similar tendency. In the figures, results for lattice size 
L/a = T/a = 6 and with the zero (BG= 0) or non-zero (BG= 1) background field 
are shown. Most eigen-values are localized near the unit circle, and a remnant of 
the GW relation is observed. Especially, note that p = is most strongly affected 
by the boundary. When switching on the background gauge field (BG= 1), some 
degeneracies are lifted, and you can see 'more' points than at zero background field 
case (BG= 0). 
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Figure 2: Spectrum of dD^ with parameters 8 = and for lattice size L/a = T/a = 6. 
The upper panels are for s = 0.0, while the lower ones are for s = 0.5. The value of BG 
means that 0: zero background gauge field, 1: non-zero background gauge field (choice 
A). All eigen- values are enclosed by a black circle whose origin is (1, 0) and radius is 1, on 
which the spectrum of the GW relation operator lie. The blue circle points which come 
from the p = sector are located far away from the circle, and they are positioned around 
a center of the circle. We observe that this sector is strongly affected by boundary effects. 



3.2 Spectrum of D^D^ 

We also investigate the spectrum of the Hermitian operator L 2 Dj <f D^ = (L/a) 2 (l + 
s) 2 aD^aD^, because we can take the continuum limit and compare the scaling 
behavior with that of the Wilson-Dirac operator, L 2 D\,D W = (L / a) 2 aD^aD^ 
|25j . The numerical results of the lowest 10 eigen- values of L 2 D^D^ with non- 
zero background field are summarized in Table [2] for s = 0.0,0.5, 9 = 0,ir/5 and 
L/a = T/a = 6, 12, 24. The scaling to the continuum limit are plotted in Figure [3] 
including those of the Wilson-Dirac and the clover action (c sw = 1) for comparison. 
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In the figure, the lower modes show a good scaling behavior, while higher modes 
are strongly affected by lattice artifacts. 



4 SF coupling to one-loop order 



4.1 Definition and results 

We compute the fermion part of the SF coupling [25] (we set L = T as usual) at 
one-loop order pn(L/a) for the massless Neuberger-Dirac operator. The one- loop 
coefficient H is given as 



1 d 

pi tl (L/a) = -— lndet£»N 



with a normalization 



ri=u=0 



1 



k = 12(L/a)^[sin(7) + sin(2 7 )], 7 = -n(a/Ly. 



(23) 



(24) 



In the actual calculation, we expand the 77 derivative and use the fact that the 
determinant is factorized to the individual spatial momentum p and color sector b, 



-Tr 



D. 



f drj 



p 6=1 L ' 



(25) 



Since Z?n is not the block tri-diagonal in the time and the spinor index, unfortunately 
we can not use the nice recurrence formula [25] which was used in the case of the 
Wilson-Dirac fermion. Therefore, we have to evaluate the inverse of -Dn directly by 
making use of a solver routine, and multiply with the 77 derivative of -Dn to take the 
trace. The trace in eq. (|25p . tr, concerns with the spinor and the time indices. 

To compute the one- loop SF coupling, we need to take the r/ derivative of -Dn. 
This can be done analytically. To this end, we have to evaluate, 



N 



f N (x) = Y J n(Y) Ck , 



(26) 



fc=0 



where the dotted defines the derivative with respect to 77. This summation can be 
evaluated by another recurrence relation besides the one needed to compose 
itself (the Clenshaw recurrence relation). The additional recurrence relation can be 
derived from the three terms recurrence formula for the Chebyshev polynomials and 
its 77 derivative formula. 



2 Formally we should define the coupling 



Pi,i(L/a) = i^lndetpDjj-DN] 



rj—b'—O 



(22) 



in order to properly define the determinant in the continuum limit, but on the lattice the above form is 
equivalent to eq. ([23| due to 75 Hermicity. 
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Figure 3: The a/L dependence of the lowest 10 eigen-values of L 2 D^D for the various 
fermion actions in the presence of the background field. The left (Right) panel is for 
8 = {6 = 7r/5). We compute the eigen- value of the Neuberger-Dirac (overlap) operator 
for s = 0.0,0.1,0.2,0.3,0.4,0.5, but here we only show s = 0.0,0.5, since the others(s = 
0.1,0.2,0.3,0.4) just interpolate in between s = 0.0 and s = 0.5. Those of the Wilson- 
Dirac action and the clover action from [25] are shown for comparison. The eigen-values 
converge to the continuum ones (Red points at a/L = 0) in Ref |25J. Within the lowest 
10 eigen-values, no level crossing occurs in the Neuberger-Dirac operator unlike the case 
of the Wilson-Dirac fermions. 
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We compute p\ : \ on the lattices of size L/a = 4, ...,48. The results are summa- 
rized in Table [3] for s parameters, s = 0.0,0.5 and 8 = 0, 7r/5. In order to estimate 
rounding off errors, we perform two methods of the approximation to the inverse 
square root as mentioned in Section \2.3\ the minimax and the Chebyshev interpola- 
tion. In the table, we show nine significant digits where both approximations agree 
with each other. Even though we have used double precision arithmetic and been 
demanding 10~ 13 precision for the inverse square root, we lose three to four digits 
in the summation step of all momentum and color sectors in eq. (|25p . 



4.2 Coefficients of Symanzik's expansion 

From the Symanzik's analysis of the cutoff dependence of Feynman diagrams on the 
lattice, one expects that the one-loop coefficient has an asymptotic expansion 

oo 

Pi,i(L/a) = Y,(a/L) n [A n + B n \n{L/a)\. (27) 
n=0 

We can reliably extract first several coefficients by making use of the method in 
Ref. [26]. 

For the usual renormalization of the coupling constant, Bq should be 2&o,i where 
&o,i is the fermion part of the one-loop coefficient of /3-function for TVf flavors QCD, 

fro = 6o,o + JVf6o,i, (28) 
fro,o = ^ (29) 



We confirmed B = 26 ,i = —0.00844343... to three or four significant digits for 
all cases (8 = 0, tt/5 and —0.6 < s < 0.6). When the tree-level O(a) improvement 
is realized, we expect that B\ = holds. We check this to 10~ 2 or 10 -4 in all 
cases. This shows that even though we have been using the approximate formula 
of the boundary coefficient in eq.©, it works well to achieve the tree- level O(a) 
improvement to the precision here. In the following analysis we set exact values 
B = -l/(12vr 2 ) and B x = 0. 

Aq gives an information about a ratio of A-parameters, and we show the obtained 
values in Table [TJ By combining the previous results from Ref. [251 EZ] ; the values 
of Aq can be obtained, and are shown in the second line (numbers with *) in each 
s in Table [1] for 6 = 0, 7r/5. We observe excellent agreements within errors for all s 
and 8 parameters we investigated. 

To achieve one-loop O(a) improvement, A% is needed to determine the coefficient 
of the fermion part of the boundary counterterm, c^ 1 ' 1 ^ [25]. The resulting values 
are shown in Table HJ No 8 dependence on the A\ is observed beyond errors. The 
absolute value of \A\\ = 0.02 — 0.03 of the Neuberger-Dirac operator is roughly 
factor two smaller than that of Wilson type fermion, \A\ \ = 0.038282(2) [25]. If one 
imposes an improvement condition [25], one finds that 



,(14) 



At/2. (31) 
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# = 


9 = tt/5 


s 


A) 


A 

Ai 


A 


A 

Ax 


u.u 


0.016944(7) 
0.016937 


-0.021(2) 


0.015562(8) 
0.015555 


-0.021(3) 


yj.o 


0.015712(6) 
0.015708 


-0.020(2) 


0.014330(5) 
0.014326 


-0.020(2) 


— n 4 


0.014754(5) 
0.014751 


-0.019(1) 


0.013373(3) 
0.013370 


-0.019(1) 


— n ^ 


0.013992(4) 
0.013990 


-0.019(1) 


0.012610(3) 
0.012609 


-0.0187(7) 




0.013385(4) 
0.013383 


-0.019(1) 


0.012003(2) 
O.OlzOOz 


-0.0187(6) 




0.012912(4) 

n m on 1 1 * 
0.01/911 


-0.019(1) 


0.011530(2) 

n ni 1 c on* 

0.011529 


-0.0188(5) 


n n 


0.012567(3) 
0.012566 


-0.0192(9) 


0.011185(1) 

n m 1 1 oc* 

0.011185 


-0.0191(4) 


n i 

U. -L 


0.012354(3) 
0.012353 


-0.0198(9) 


0.010972(1) 
0.010972 


-0.0197(4) 


n 9 


0.012287(3) 
0.012285 


-0.0207(9) 


0.010905(2) 
0.010904 


-0.0206(4) 


n 3 

u.o 


0.012390(3) 

n m 9^88* 

U.Ul ZjOOO 


-0.0222(8) 


0.011008(2) 

n m 1 007* 

U.Ul 1UU I 


-0.0220(4) 


0.4 


0.012704(3) 
0.012703* 


-0.024(1) 


0.011322(2) 
0.011321* 


-0.0241(6) 


0.5 


0.013293(7) 
0.013292* 


-0.028(2) 


0.011912(3) 
0.011911* 


-0.0275(9) 


0.6 


0.01426(2) 
0.014259* 


-0.032(7) 


0.012880(3) 
0.012878* 


-0.033(1) 



Table 1: The coefficients of asymptotic expansion for —0.6 < s < 0.6 for 9 = 0, 7r/5. 
The value of Aq with * in the last line in each s parameter block are the values from the 
previous calculations [25l [27]. The error for those values should be on the last digit. 

For future reference, we provide interpolation formula of the c^ 1 ' 1 ^ as a polynomial 
of s for 9 = 7r/5 where errors are smaller than 9 = case 

cj 1 ' 1 ) = -0.00958 - 0.00206s - 0.00484s 2 - 0.00748s 3 - 0.01730s 4 , (32) 
for -0.6 < s < 0.6. 
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5 Lattice artifacts of the step scaling function 
to one-loop order 

In this section, we investigate lattice artifacts of the step scaling function (SSF) 
|28| a(2,u), which describes the evolution of the running coupling g 2 (L) = u under 
changes of scale L by a factor 2, 

a(2,u) =g\2L), u = g 2 {L). (33) 

The lattice version of the step scaling function is denoted by Y,(2,u,a/L). 

Perturbative estimate of the lattice artifacts of the step scaling function can be 
studied by expanding a relative deviation 

S(u,a L)= v 11 v ; =6i (a L)u + 0{u 2 . 34 
a(2,u) 

The one- loop deviation, 5i(s,a/L), is decomposed into pure gauge and fermion part 

' ==! ' Si(a/L) = S Xja (a/L) + N { 5 1A (a/L). (35) 

We are currently only interested in the fermion part. The fermion part of the one- 
loop deviation 5i t \(a/L) in terms of the one-loop coefficient of the SF coupling p\ t \ 
is given by 

S 1A (a/L) = Pl ,i(2L/a) - pi,i(L/a) - 2b ,i ln(2). (36) 

Depending on the value of the boundary counter term c^ 1,1 ^ , we denote with s\ °\ (a / L) 
as the tree level 0(a) improved version with c t (1,1) = 0, and 8f[{a/L) the one-loop 

0(a) improved one for c^ 1 ' 1 ^ = A\/2. 

We show numerical results of the one-loop deviation in Table d] and plots in 
Figure [U where we include those of the Wilson-Dirac and the clover action for 
comparison [25]. In the case of the clover action, cj; 1 ' 1 ^ is set to be the proper value 
to achieve one- loop 0(a) improvement, and for the Wilson fermion it is set that 
c^ 1 ' 1 ^ = 0. We observe that the lattice artifacts for the Neuberger-Dirac operator 
are comparable to those of the clover action. As in the case of the clover action, 
the Neuberger-Dirac operator has less lattice artifacts for the case of 9 = ir/5 than 

6 = 0. 



6 Conclusion and outlook 

In this paper, we have explored the free Neuberger-Dirac operator on the SF. We 
investigated the spectrum of the operator, and then we confirmed that the spectrum 
of .Dn is enclosed by the unit circle and the spectrum of D^D^ has the expected 
scaling behavior (1/L 2 ) and the correct continuum limit. We also performed the one- 
loop computation of the SF coupling by making use of the operator. We confirmed 
the universality, and the fermion part of the 0(a) boundary counterterm at one- 
loop order, c^ 1,1 ^ is determined. The formula in eq. (f32j) . c^ 1 ' 1 ^ as a function of s, 
might be useful for future simulations. By making use of the one- loop results, we 
estimated the lattice artifacts of the SSF. It turns out that the size of the lattice 
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6=0 



0.005 
0.004 
0.003 
0.002 
0.001 


-0.001 
-0.002 
-0.003 
-0.004 
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- Wilson — -r^ _ 
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0.15 

a/L 

e=7t/5 



0.2 0.25 





0.07 



0.03 0.04 

(a/L) 2 



0.07 



Figure 4: We show the relative deviation with the various actions for tree level 0(a) 
improvement, 5^1 (Left), and one-loop 0(a) improvement, 5^1 (Right), as a function of 
a/L and (a/L) 2 respectively. Upper part is for 8 = 0, and lower is for 6 = it/ 5. For 
comparison, those of the Wilson-Dirac fermion with c^ 1 ' 1 ^ = and the clover fermion with 
c\ = 0.019141 [25] are included in the plot of o\ { and o\ [ respectively. 



artifacts for the Neuberger-Dirac operator is almost the same as that of the clover 
action. Thus, we may expect small lattice artifacts for the non-perturbative SSF of 
the Neuberger-Dirac operator, as in the case of the clover action [llj . In Appendix 
IB"1 we demonstrate that the choice of the boundary coefficient c = 1 + s given in 
|20| at tree level is almost optimal. This formula may be precise enough for actual 
simulations. 

We exclusively considered the massless case. Although the massive case can be 
explored, we leave it as a future task. Comparison of scaling behavior with the 
other formalism [19j and a consistency check is also interesting. Before starting 
non-perturbative computations, we have to compute some improvement coefficients 
to one or two loop order. Furthermore perturbative calculations of the renormal- 
ization factors are in a to-do list. By combining techniques in Ref. [29], a two-loop 
calculation of the SF coupling including the Neuberger-Dirac operator as a fermion 
part may also be feasible. 

Apart from the one-loop computations, next target would be a computation of 
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the renormalization constant of the flavor singlet scalar density Z$ non-perturbatively, 
since the bare quark condensate in two flavor QCD was already computed by JLQCD 
[ipl . Due to the chiral symmetry, Zg is identical to the renormalization constant 
for the non-singlet flavor pseudo scalar density, Zp = Z$ [301 [31]. Actually, the 
non-perturbative renormalization group running of Zp is already known in Ref. [12] 
for the SF scheme. A missing piece to obtain the renormalization group invariant 
quark condensate is a low energy matching factor, Zp(go,[i = 1/L max ) in the SF 
scheme for the overlap fermion. This is an urgent and possibly doable task in the 
near future. 
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A Free Wilson-Dirac operator on the SF in 
time- momentum representation with non-zero 
background field 

In the presence of the background gauge field [6] , the free part of the Wilson fermion 
action in time-momentum space has a form |25j 

S ™ = 7~3 X X ^(~P> x o)(£>w(p; x ,y ) + 5 XOt y m )ip(j>, yo), (37) 
p xo,yo 

with the boundary conditions 

P+^( P ,0) = 0, P_V(P,T)=0, (38) 
^(p,0)P_ = 0, ?P( P ,T)P + = 0. (39) 

The massless part of the Wilson-Dirac operator on the SF is given as 

aD^(p;x ,yo) = {-PS Xo+am + h b (p;x )5 XO: y - P + 5 Xo - a>yo }S bc , (40) 

where indices b, c refer to color and h b (p; xq) is given by 

3 3 

h\p; Xo ) = l + ~J2 q b k {x,f + i £ ql(x ) lk , (41) 

k=l k=l 

3 They performed a non-perturbative renormalization by the RI/MOM scheme. In this paper here, we 
are talking about the non-perturbative renormalization by making use of SF scheme. 
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with 



q b (x ) = sm(aq b (x )), (42) 

q b (x ) = 2sm(aq b (x )/2), (43) 

and 

q b k (x ) = w b x + r b k , (44) 

w b = (ti-^/L 2 , (45) 

r b k = Pk + h/L. (46) 

The spatial component of the momentum p is given by 

p k = (2vrn fc + 9)/ L, (47) 

for n k = 0, ■ ■ ■ , L/a — 1. The boundary phases (f) b an d 4>' b are given by 

, / 47T . , . 

01 = v - 3, 0i = -0i - -3-' ( 48 ) 

1 27T 

02 = r/(-- + ^), ^ = -0 3 + _, (49) 

1 7T 2TT 

03 = r/(-- - v) + -, ^ = - ( j )2 + —. (50) 

Switching off the phases <p b = (p' b = correspond to the zero background field. 
One can consider the effect of the non-zero background field a shift for the spatial 
momentum. The 9 also affect as a constant shift for the momentum in a global 
manner, on the other hands, the background field provides a time dependent (local) 
shift. 

In a 4(T/a — 1) x 4(T/a — 1) matrix expression, the Wilson-Dirac operator is 
represented by 

aD b w (p) (51) 

h b (p;a) -P- 

-P+ h b (p;2a) -P_ ••• ••• 

-P+ h b {p;3a) -P_ ••• : 

; -p+ ••. ■•. ■•. ; 

: '-. '■■ "•• -P_ 

••• ••• -P+ h b {p-T-2a) -P_ 

-P+ h b (p;T-a) 



B Determination of boundary coefficient c at 
the tree level 

The boundary coefficient c in the kernel of the overlap operator on the SF is ex- 
panded in terms of the coupling constant , 

c(gl) = c^+c^gl + 0(gt). (52) 
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In this appendix, we determine the tree coefficient which depends on s param- 
eter. Ref. [20j gives the formula in eq.([7|), and we will examine it carefully here. 

We consider the SF with T = 2L and 9 = in the presence of the non-zero 
background field. The massless Neuberger-Dirac operator is assumed also in this 
appendix. The basic correlation functions [32] we use are given by 



f A (x ) = -a 6 ^i(^(x)C(y)75^r a C(z)}, (53) 
f P (x ) = -a 6 ^i(P a (x)C(y)75^r a C(z)), (54) 
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where boundary fields [20j are given by 



C(x) = U(x,0)\ Xo=0 P^(x)\ Xo=a , (55) 
C(x) = ^(x) U 0=o ^+i7(x, 0)- 1 U o=0 . (56) 



At the tree level, /a(^o) fp( x o) are given as 



fPHxo) = £/S(zo), (57) 



a=l 



f?l( x o) = -tv[P + TP^S a (p = 0;a,x )TS a (p = 0;x ,a)}, (58) 



with T = 7075,75 for and f p °^ respectiveljQ, and 5 a (p;xo,yo) is a free propa- 
gator in time-momentum space. The trace in the above equation are over the Dirac 
spinor indices only. The a refers to color. We chose a ratio f^\xo)/fp\xo) at 
a middle point xo = T/2 as an observable to impose the improvement condition. 
We compute the quantity ff ] {T/2)/ f { p 0) (T/2) for lattice size L/a = 8, 10, ...,64 and 
—0.6 < s < 0.6 and some range of . We search around the target point with 
width Ac(°) = 0.0005. 

We extract the order a coefficient A\ from the Symanzik's expansion for the 
ratio 

f (0) Irp ln\ °° 

%™ = J>/L)M n . (59) 
fP(T/2) fo 

We estimate the error of A± by making use of the method in |26j. 

We determine c*^ such that Ai(c*^) = for the range —0.6 < s < 0.6 (im- 
provement condition). In Figure [SJ we plot c*(°'(s) as a function of s. By fitting 
the data points with a functional form 

c *(°)(s) = 1 + kis + k 2 s 2 + k 3 s 3 + hs* + hs 5 , (60) 

we obtain 

h = 1.0002, k 2 = -0.1279, k 3 = -0.0374, k 4 = 0.1616, k 5 = 0.0818. (61) 



One has to use improved operators, and this is equivalent to replace the propagator D^ 1 — > (1 — 
aD /2)D~ 1 in the case of massless. The second term turns out to be a contact term, and drops when one 
considers a correlation function whose insertion points are separated like here. 
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This curve is also shown in Figure [5j For larger s, a discrepancy between the above 
formula and eq.© can be seen, and their difference is maximally 10% in the range 
—0.6 < s < 0.6. As a consistency check, by making use of the value of c in eq. (|6UI6ip . 
we compute the ratio in eq . (|59f) with 9 = vr/5, and then A\ = is confirmed up to 
10 -4 in the range of s. 




-0.6 -0.4 -0.2 0.2 0.4 0.6 



Figure 5: The green dashed line represent c*'°'(s) in eq. (l60|) . although * is not shown in 
the legend. The line of c(s) = 1 + s which is given in Ref. [20] is also shown as solid red. 
The error bar of the points are too small to see in this scale. 

In order to measure an effect of the difference on a physical quantity, we compare 
the onedoop coefficient px^ with different values of c from the different formulae 
of c (eq.© and eq. ([60]) ) at s = 0.5. It turns out that a difference in the px t i is 
less than one percent on the lattice size L/a = 4, ...,48. Furthermore, the resulting 
Symanzik's coefficients of pn in eq. (|27|) do not change within errors. Therefore, we 
concluded that, to the precision in our calculation, the formula c = 1 + s is accurate 
enough to achieve the tree level O(a) improvement. 

C Tables of numerical results 
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e 


= 












L/a 


= 6 


L/a 


= 12 


L/a 


= 24 






n 


s = 0.0 


s = 0.5 


s = 0.0 


s = 0.5 


s = 0.0 


s = 0.5 


b 


d 


1 


2.707535 


2.544076 


2.353314 


2.293247 


2.228076 


2.202150 


2 


2 


2 


5.443956 


5.303784 


5.178180 


5.079915 


5.002836 


4.946708 


2 


2 


3 


9.738741 


9.383681 


8.485551 


8.277036 


7.995129 
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3 


2 


4 
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10.185393 


1 


2 


5 
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3 


2 


6 
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20.658775 


1 


2 


7 


25.715823 


26.315166 


24.501908 


24.335325 


23.730591 


23.572580 
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2 


8 


26.985530 


27.522077 


25.804837 


25.514692 


24.860285 


24.610036 


2 


2 


9 
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34.128759 
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29.464037 


28.196362 


28.035155 


1 


6 


10 


32.231451 


34.580806 


30.746251 


30.712437 


29.585099 


29.393814 


3 


6 








9 = 


= tt/5 












L/a 


= 6 


L/a 


= 12 


L/a 


= 24 






n 


s = 0.0 


s = 0.5 


s = 0.0 


s = 0.5 


s = 0.0 


s = 0.5 


b 


d 


1 


6.011926 


5.711491 


5.218218 


5.083175 


4.924877 


4.861719 


2 


2 


2 
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6.262644 
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5.082355 


1 


2 


3 
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2 


2 


4 
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14.477736 
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13.748289 
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13.413035 


1 


2 


5 
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3 


2 


6 


22.198468 
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3 
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7 
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8 
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28.332464 
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2 


2 


9 
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28.909787 


28.721464 


1 


6 


10 


32.344003 


34.955020 


30.320555 


30.299555 


28.987432 


28.807336 


3 


6 



Table 2: The lowest 10 eigen-values of the Hermitian operator L 2 D\^D^ for s = 0.0, 0.5. 
Upper (Lower) panel is for 9 = (9 — 7r/5). The b represents the color sector, and the d 
is for degeneracy for one flavor. 
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Table 3: The one loop coefficient pi : i(L/a) for s = 0.0, 0.5 with 9 = 0, n/5. The last digits 
may be affected by rounding off errors. 
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00079 





00154 





00262 


-0 


00264 


-0 


00046 


0.00034 


0.00087 


0.00058 


9 





00082 





00134 





00212 


-0 


00270 


-0 


00029 


0.00027 


0.00057 


0.00041 


10 





00082 





00116 





00176 


-0 


00268 


-0 


00018 


0.00020 


0.00036 


0.00030 


11 





00079 





00102 





00151 


-0 


00263 


-0 


00012 


0.00015 


0.00024 


0.00023 


12 





00075 





00091 





00132 


-0 


00256 


-0 


00008 


0.00011 


0.00016 


0.00018 


13 





00071 





00082 





00118 


-0 


00249 


-0 


00006 


0.00008 


0.00011 


0.00015 


14 





00067 





00075 





00108 


-0 


00241 


-0 


00004 


0.00006 


0.00008 


0.00012 


15 





00063 





00069 





00099 


-0 


00234 


-0 


00004 


0.00005 


0.00006 


0.00011 


16 





00059 





00064 





00092 


-0 


00226 


-0 


00003 


0.00004 


0.00005 


0.00009 


17 





00056 





00059 





00086 


-0 


00219 


-0 


00003 


0.00003 


0.00004 


0.00008 


18 





00053 





00056 





00081 


-0 


00213 


-0 


00003 


0.00002 


0.00003 


0.00007 


19 





00050 





00052 





00076 


-0 


00206 


-0 


00003 


0.00002 


0.00003 


0.00006 


20 





00048 





00049 





00072 


-0 


00200 


-0 


00002 


0.00001 


0.00002 


0.00006 


21 





00045 





00047 





00069 


-0 


00195 


-0 


00002 


0.00001 


0.00002 


0.00005 


22 





00043 





00045 





00066 


-0 


00189 


-0 


00002 


0.00001 


0.00002 


0.00005 


23 





00041 





00042 





00063 


-0 


00184 


-0 


00002 


0.00001 


0.00002 


0.00004 


24 





00040 





00041 





00060 


-0 


00179 


-0 


00002 


0.00001 


0.00002 


0.00004 


9 = tt/5 




"1,1 


°1 1 


L/a 


s 


= -0.5 


s 


= 0.0 


s 


= 0.5 


Wilson 


s - 


= -0.5 


s = 0.0 


s = 0.5 


Clover 


4 





00072 





00241 





00475 


-0 


00273 


-0 


00178 


0.00002 


0.00132 


0.00009 


5 





00058 





00196 





00360 


-0 


00330 


-0 


00142 


0.00005 


0.00085 


-0.00005 


6 





00067 





00148 





00254 


-0 


00346 


-0 


00100 


-0.00011 


0.00025 


-0.00010 


7 





00071 





00115 





00191 


-0 


00344 


-0 


00071 


-0.00021 


-0.00005 


-0.00010 


8 





00071 





00095 





00155 


-0 


00334 


-0 


00054 


-0.00024 


-0.00017 


-0.00008 


9 





00067 





00083 





00135 


-0 


00322 


-0 


00044 


-0.00024 


-0.00018 


-0.00007 


10 





00062 





00074 





00121 


-0 


00309 


-0 


00038 


-0.00021 


-0.00017 


-0.00006 


11 





00058 





00068 





00111 


-0 


00296 


-0 


00033 


-0.00019 


-0.00014 


-0.00005 


12 





00054 





00063 





00102 


-0 


00284 


-0 


00029 


-0.00017 


-0.00012 


-0.00004 


13 





00051 





00058 





00096 


-0 


00273 


-0 


00026 


-0.00015 


-0.00010 


-0.00003 


14 





00048 





00055 





00090 


-0 


00263 


-0 


00023 


-0.00013 


-0.00009 


-0.00003 


15 





00046 





00052 





00084 


-0 


00253 


-0 


00021 


-0.00012 


-0.00007 


-0.00003 


16 





00044 





00049 





00079 


-0 


00244 


-0 


00019 


-0.00011 


-0.00007 


-0.00002 


17 





00042 





00047 





00075 


-0 


00235 


-0 


00017 


-0.00010 


-0.00006 


-0.00002 


18 





00040 





00044 





00071 


-0 


00227 


-0 


00015 


-0.00009 


-0.00005 


-0.00002 


19 





00039 





00042 





00068 


-0 


00220 


-0 


00014 


-0.00008 


-0.00005 


-0.00002 


20 





00037 





00040 





00064 


-0 


00213 


-0 


00013 


-0.00007 


-0.00004 


-0.00002 


21 





00036 





00039 





00062 


-0 


00206 


-0 


00012 


-0.00007 


-0.00004 


-0.00001 


22 





00035 





00037 





00059 


-0 


00200 


-0 


00011 


-0.00006 


-0.00004 


-0.00001 


23 





00033 





00036 





00056 


-0 


00194 


-0 


00010 


-0.00006 


-0.00003 


-0.00001 


24 





00032 





00035 





00054 


-0 


00189 


-0 


00009 


-0.00005 


-0.00003 


-0.00001 



Table 4: The relative deviation 5[ { for s = —0.5, 0.0, 0.5 for 6 = O(upper) and 
6 = 7r/5(lower). Results for the Wilson and the clover action are also included 
for comparison. 
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